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ABSTRACT 

We present Suzaku observations of the Centaurus cluster out to 0.95r2oo, taken along a strip to 
the north west. We have also used congruent Chandra observations of the outskirts to resolve 
point sources down to a threshold flux around 7 times lower than that achievable with just 
Suzaku data, considerably reducing the systematic uncertainties in the cosmic X-ray back- 
ed ■ ground emission in the outskirts. We find that the temperature decreases by a factor of 2 
O I from the peak temperature to the outskirts. The entropy profile demonstrates a central excess 

(within 0.5r2oo) over the baseline entropy profile predicted by simulations of purely gravi- 
tational hierarchical structure formation. In the outskirts the entropy profile is in reasonable 
agreement with the baseline entropy profile from Voit et al., but lies slightly below it. We 
find that the pressure profile agrees with the universal pressure profile of Arnaud et al. but 
lies slightly above it in the outskirts. The excess pressure and decrement in entropy in the 
outskirts appear to be the result of an excess in the measured gas density, possible due to gas 
clumping biasing the density measurements high. The gas mass fraction rises and reaches the 
mean cosmic baryon fraction at the largest radius studied. The clumping corrected gas mass 
fraction agrees with the expected hot gas fraction and with the simulations of Young et al. We 
■ further the analysis of Walker et al. which studied the shapes of the entropy profiles of the 

clusters so far explored in the outskirts with Suzaku. When scaled by the self similar entropy 
the Suzaku entropy profiles demonstrate a central excess over the baseline entropy profile, 
and are consistent with it at around t^qq. However outside r5oo the entropy profiles tend to lie 
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\ below the baseline entropy profile. 
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1 INTRODUCTION temperature and density profiles for data from within rsoo ~ §^200, 

. (for example Vikhli nin et al] 2006), meaning that only ~ 30 percent 
The low and stable particle background of Suzaku has allowed sig- . ^ , , , ... , , _, . . , . 

, ,.f r , , r , ■ , oi the cluster volume had actually been explored. This is subject to 

niricant progress to be made in the study or the low surface bright- . , . .... . . 

large errors and priors, as there is no robust prior expectation tor 
ness X-ray emission from galaxy cluster outskirts, providing ex- c . , , „, . . , 

J r ° the form of the temperature and density profiles in the outskirts, 

citing and unique observations of these previously unexplored re- 
gions. In the outskirts, gas is continuing to accrete onto clusters, 



Accurate cluster mass determinations are important for using 

clusters as probes of cosmological parameters using the mass func- 
allowing us to see the formation process in action. We can gain an . . u... , ,. . — t , \\~ nn A ITT — — 711™, A, T 

6 „ * 6 Hon method ( Vikhlmin et al. 2006, Mantz et al. 2010). Improving 



our understanding of the gas mass fraction of clusters at large ra- 
dius allows us to better understand the assumptions behin d usi ng 



understanding of how close the ICM is to hydrostatic equilibrium 
and spherical symmetry in the outskirts, which are the fundamental 

assumptions made when deriving masses of galaxy clusters from , ,, . , -,,„,.-, 

rn the gas mass traction as a cosmological probe (Allen et al. 2008) 
X-ray observations. Estimations of the cluster mass within r2od_l ^, i 11 1 ■ 1 , r 1 * ■ ■ ^ • 1 o , 

have previously involved extrapolating outwards analytic best fit 



The breakthrough in the study of cluster outskirts with Suzaku 
was made for clusters whose size and redshift allowed the virial 
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1 f200 is the radius within which the mean density of the cluster is 200 cally used to represent the virial radius. The mass enclosed within r200 is 

times the critical density required for a flat universe, p c , and is typi- M200=4/37rp c r2QQ. 
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radius to be reached using a small number of pointings, in that 
the angular extent of the cluster made the distance between the 
core and r2oo around 20-30 arcmins. These clus ters are PKS 
0745-191 dGeorge et al.l 120091: IWaJkeretal.ll2012al) Abell 2204 
Reiprich et alJ|2009h . Abell 1795 dBautz et al.ll2009h Abe ll 1413 
ashino et alj|2010h. Abell 1689 dKawaharada et alj|20ich . Abell 
2142 jAkamatsu et all l20lTh . th e fossil group R X Jl 159+5531 
( Hum phrey et alj2012D . Hydra A dSato et al.l2012l) and Abell 2029 
(Walk eret al.l2012l) . However the large point spread function (PSF) 
of Suzaku (HPD=2') limits the size of annuli used in spatially re- 
solved spectral analysis, reducing the resolution of the temperature 
and entropy profiles. 

To get around the large PSF and improve the spatial resolution 
of the profiles of the ICM properties, we need to observe nearby , 
X-r ay bright clusters, as was achieved in lSimionescu et all d201 ll) 
and lSimionescu et "all d2012h for the Perseus cluster. The Centau- 
ries cluster's low redshift (z=0.0109) gives an angular extent of 13 
kpc/arcmin, meaning more spatially resolved temperature and en- 
tropy profiles can be obtained. Centaurus also lies further from the 
galactic plane than Perseus giving it a lower absorbing column and 
reducing the fluctuations in column density across the face of clus- 
ter. 

At present, Suzaku has mainly studied high mass clusters to 
the virial radius [the fossil group RX J 1 159+5531, Hydra A (3keV) 
and Virgo (2.3keV) are the only clusters below 4 keV to be studied 
to r2oo]. Centaurus' low average temperature (we find fcT(0.1 < 
r < 0.5r 200 ) = 3.0±£;i keV) means that we are probing the low 
mass range of clusters, and so its study is of importance as it will 
allow us to gain an understanding of the properties of clusters of 
different masses in the outskirts. It is important to explore the full 
mass range of clusters to investigate M — T,M — L and Yx — M 
scaling relations t hroughout the complete mass range. For instance, 
lArnaud et alJfcOOST) found that when clusters below kT (0.1 < r < 
0.5r2oo) = 3.5 keV are included in the M-T relation, it steepens 
compared to the self-similar relation M cc T , becoming M cc 

rpl.7 

Of particular interest is the entropy profile (K = kT/rij^ 3 ). 
Nume rical simulations of pure gravitational collapse dVoit et al.l 
2005) predict that the entropy should increase as a powerlaw with 
radius (K oc r ) and also provides a prediction for the normalisa- 
tion of this baseline entropy profile. Deviations from this baseline 
entropy profile must originate from non-gravitational processes, 
and so comparing the shape and normalisation of the observed en- 
tropy profile with the baseline entropy profile provides an insight 
into the physical processes occurring i n the ICM. 

We have already shown in Wal ker et al.l (2012b) that the en- 
tropy profiles for clusters explored with Suzaku out to r2oo have 
the same shape, flattening away from a powerlaw increase above 
~0.5r5oo- The same entropy profile shape was also found f or the 
XMM-Newton study of the Virg o cluster dUrban et al.ll20lH) . and 
the Chandra studies of A2204 dSanders et alj 120091) and A 1835 
dBonamente et al.ll2013l ). However, more insight can be obtained 
by scaling the entropy profiles by the self similar entropy at rsoo, 
which allows both that shape and normalisation to be compared to 

the baseline profile. This is explored in section|6] 

A similar entropy analysis has been performed in IPratt et all 

d2010l) for the entropy profiles of the REXCESS clusters (XMM- 
Newton observations of 31 cluste rs extending to a t least riooo, and 
up to rsoo for 13 clusters) and in lSun et"al] d2009h for a sample of 
groups. In both cases the entropy profiles are flatter than the K oc 
r 1,1 powerlaw prediction at rsoo- However the entropy is found to 
exceed the baseline level inside rsoo, such that the flattening only 



acts to bring the entropy level into agreement with the baseline level 
at rsoo, without ever requiring it to go below the baseline entropy 
profile. This entropy excess is possibly the result of a combination 
of extra heating and mixing of the ICM caused by mergers. The 
central entropy excess is mass dependent, with less massive clusters 
showing a greater excess. 

Also of interest is the pressure profile. lArnaud et al.l (|2010|) 
used the REXCESS sample of clusters to show that the pressure 
profiles of clusters tend to have the same shape and normalisa- 
tion when scaled by the characteristic pressure, P500, and have 
used predictions of simulations to extend the proposed universal 
pressure profile out to and beyon d r2oo • Observations with Planck 
( Planck Collaboratio n et al .120121) using the Sunyaev-Zeldovich ef- 
fect to measure the pressure profiles of a stacked sample of 62 clus- 
ters have agreed with the form of this universal pressure profile in 
the outskirts of clusters, but has found slightly higher pressures than 
simulations predict outside rsoo. 

Cluster outskirts observations are also important in constrain- 
ing the level of gas clumping o ccurring the ICM, whic h is a predic- 
tion of numerical simulations dRoncarelli et al.ll2006l) . Gas clump- 
ing causes the gas density to be overestimated if the ICM is as- 
sumed to be unif orm, resulting in underestimates of the entropy 
dNagai & Laii201 lh . and overestimates of the gas pressure and gas 
mass fraction, as has been observed for Perseus dSimionescu et al.l 
l201lh and PKS 0745-191 I Wal ker et alj|2012al) . Simulations also 
predict that the contribution of non-thermal pressure suppor t in the 
ICM increases to around 20 percent at r2oo (Lau et al. 2009), caus- 
ing hydrostatic equilibrium masses which use only the gas pressure 
to underestimate the total cluster mass. 

We can also study the gas mass fraction profile, which de- 
pends on cluster m ass, and which will allow feedback models 
dYoung et al JI201 lh to be tested. In less massive clusters it is ex- 
pected that feedback processes are more able to redistribute gas to 
larger radius, and these observations provide an excellent way of 
testing these models. 

Here we present Suzaku observations of the Centuarus cluster, 
taken along a strip from the core to the outskirts in the north west 
direction, reaching out to r2oo and beyond. The north west direction 
was chosen bec ause it avoids the c old front to the west (studied 
with Chandra in Fabian et al. l2005h . and also avoids the merging 
activity reported in IChurazov et al.1 d 19991) to the south east. We 
therefore expect this direction to be relaxed, allowing a hydrostatic 
mass analysis to be performed and a determination of the gas mass 
fraction. 

We use a standard ACDM cosmology with Ho = 70 km s _1 
Mpc~\ tl M = 0.3, fi A =0.7. All errors unless otherwise stated are 
at the 1 a level. 

All spectral fits were performed in XSPEC 12.7.0u using the 
extended C-statistic. 



2 OBSERVATIONS AND DATA REDUCTION 

Six observations were taken in a strip to the north west, avoiding 
known cold fronts and bright point sources. The roll angle was cho- 
sen so that pointing diagonal points towards the core, which has 
been found to minimise the effects of stray light (as shown in Figs. 
6.18 and 6. 19 of the Suzaku Technical DescriptiorQ) due to the de- 
sign of the X-ray telescopes (which were built as four sectors). The 



2 heasarc.nasa.gov/docs/suzaku/prop_tools/suzaku_td/suzaku_td.html 
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Figure 1. Le/hMosaic of Suzaku pointings in the 0.7-7.0 keV band. The annuli used in spectral extraction are shown in green and the ring radii are at 
5, 10,15, 20,25, 30,35, 45, 50, 60,70, 80,and 90 arcmins. The value of r'200 calculated later in the mass analysis is shown by the red line. The cyan sector shows 
the region which is free from stray light due to the pointings being orientated with the chip diagonal towards the core. Point sources which are resolved in the 
Suzaku images have been removed from the image (and from the spectral analysis). Centre: Mosaic of Chandra ACIS-I observations (0.5-7.0 keV band) of 
the regions more the 50 arcmins from the cluster centre. Point sources identified with WAVDETECT are circled. The spectral extraction annuli are the same as 
shown in the left panel. Right: Mosaic of ROSAT PSPC pointings with the Suzaku pointing locations overlaid in red. The red circle shows the value of r2oo 
we calculate later. The offset pointing is of NGC 4507 and is used as a background pointing in the ROSAT analysis we perform. 
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stray light free region is a sector of angle 12.8 degrees when the 
chip diagonal points towards the cluster core (see Fig. 6.19 of the 
Suzaku Technical Description), and this is shown by the cyan sector 
in Fig.Q] which shows that in the outskirts the pointings are largely 
free of stray light. We use an archival Suzaku core observation to 
complete the mosaic image. The observations used are shown in 
tableQ] and the mosaic image of these pointings is shown in Fig.[T] 
left. 

Four lOks Chandra ACIS-I observations were obtained (PI: 
S.A. Walker) of the regions covered by the outermost Suzaku point- 
ings (outside 50 arcmins from the core), to better resolve point 
sources and reduce the systematic errors of the uncertainty of the 
CXB level on the measurements. The importance of using Chandra 
to resolve point sources in Suzaku a nalyses of cluster outskirts is 
also described in lMilleretal.N2012h . The observations are shown 
in Fig. Q] (centre) and the details are shown in table Q] 

In addition we use 7 archival ROSAT PSPC observations 
shown in table Q] and mosaicked in Fig. Q] (right), which are used 
later in determining the expected spatial variations of the soft galac- 
tic background components, and in calculating a density profile to 
check the Suzaku results. 



Figure 2. Checking for SWCX contamination during the periods of obser- 
vation using the proton flux obtained with the WIND spacecraft. SWCX 
is important when the proton flux is greater then 4 X 10 s cm~ 2 s~ , Verti- 
cal dashed lines show the start and end of each observation. The pointings 
where taken in order of increasing radius from the centre. The only point- 
ing for which the proton flux exceeds 4x 10 s cm -2 s — 1 is the central most 
pointing (806084010). The fact that the excess is only slight, combined with 
the high brightness of the cluster in this region near the core, makes the ef- 
fects of possible SWCX contamination negligible. 



2.1 Suzaku data reduction 

The Suzaku data w ere reduced using the method described in 
IWalker et all d2012h using HEAsoft version 6.12 and the CALDB 
released on 2012 October 5, and we used the latest contamination 
layer calibration. The observations were reprocessed using the ftool 
AEPIPELINE, which performs the standard cleaning described in ta- 
bles 6.1 and 6.2 of the Suzaku ABC guidJfL and in addition we 
used the COR>6 condition. Spectra were then extracted from the 



3 http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/abc/ 
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Table 1. Observational parameters of the pointings 



Instrument 


Obs. ID 


Position 


Total exposure 


RA 


Dec (J2000) 


Date 


N H 1 10 20 cm" 2 








per detector (ks) 








Kalberla et al. (20051 


Suzaku XIS 


800014010 


Centre 


28.2 


192.2012 


-41.3132 


2005-12-27 


8.30 




806084010 


1 


10.9 


192.0509 


-41.0558 


2012-01-16 


8.43 




806085010 


2 


13.4 


191.8847 


-40.8272 


2012-01-16 


8.41 




806086010 


3 


10.0 


191.7105 


-40.5689 


2012-01-16 


8.67 




806087010 


4 


10.1 


191.5735 


-40.3305 


2012-01-17 


8.05 




806088010 


5 


13.2 


191.4290 


-40.0920 


2012-01-17 


8.14 




806089010 


6 


11.1 


191.2847 


-39.8559 


2012-01-17 


7.76 


Chandra ACIS-I 


15182 




10 


191.675 


-40.448 


2012-11-26 






15183 




10 


191.488 


-40.18 


2012-11-26 






15184 




10 


191.324 


-39.91 


2012-11-25 






15185 




10 


191.1475 


-39.638 


2012-11-25 




ROSAT PSPC 


rp800192n()0 


Centre 


7.793 


192.2000 


-41.31000 








rp800607n00 


Centre 


6.787 


192.2000 


-41.31000 








rp800323n00 


North 


20.658 


192.5500 


-40.52000 








rp800321n00 


East 


1.584 


193.2600 


-41.57000 








rp800322n00 


West 


16.908 


191.1500 


-41.06000 








rp800324n00 


South 


3.782 


191.8600 


-42.11000 








rp701518n00 


NGC 4507 


5.882 


188.9000 


-39.91000 







XIS0, XIS 1 and XIS3 detectors using the annulus regions shown in 
Fig- E] The calibration regions at the corners of the detectors were 
removed, and the regions at the edges of the detectors where the 
effective area is low were removed. For the XIS0 detector, the de- 
fective region (from a micrometeorite hit) was removed. Outside 
r — 50 arcmins, the small regions of the Suzaku data which did 
not overlap with the Chandra data were excluded from the spectral 
extraction. 

For each spectrum the non X-ray background (NXB) spectrum 
was produced using the ftool XISNXBGEN which uses a database of 
night earth observations, and we used the standard 300 day interval 
(between 150 days before and 150 days after the observation) to 
construct the NXB spectra. The latest calibration files were used, 
which account for the modification of the NXB level of the XIS 1 
detector due to the increase in the charge injection level on 2011 
June 1. We produced two ARFs for each spectrum using XISSI- 
MARFGEN due to the different spatial variations of the background 
and cluster emission. The first assumes a uniform source and is 
used for the background model when performing spectral fitting 
in xspec. The second used a background subtracted image of the 
cluster, and is used for the model of the cluster emission. 

The light curves were checked to ensure that no flaring oc- 
curred during the remaining cleaned observations. We then checked 
that the observations were not contaminated by solar wind charge 
exchange emission (SWCX) by investigating the proton flux from 
the WIND spacecraft's SWE (Solar Wind Experiment) instru- 
men t as shown i n Fig. [2] It has been found dFuiimoto et alj2007l 
Yoshin o et al.l2009h that Suzaku spectra do not show strong SWCX 
signatures when the proton flux is below 4x10 s cm -2 s _1 . As 
shown in Fig. [2] the proton flux only slightly exceeds 4x 10 s cm -2 
s _1 during the first observation (806084010), which is the closest 
to the cluster core, while the other observations are not affected. 
The high X-ray brightness of the ICM emission from this pointing, 
together with the fact that this is only a slight excess in the proton 



4 http://web.mit.edU/afs/athena/org/s/space/www/wind.html 



flux (the proton flux can increase by an order of magnitude above 
4x 10 8 cm -2 s -1 during a strong flare), makes the effect of SWCX 
negligible. 

The mosaicked Suzaku image shown in F ig. [TJ was obtained 
following the method described in lBautz et al.l l|2009), adding to- 
gether the images from the front illuminated detectors (XIS0, 
XIS3). Point sources which are bright enough to be resolved by 
the Suzaku pointings have been removed from this image and from 
the spectral extraction using 2.5 arcmin circular exclusion regions. 

2.2 Chandra data reduction 

The ACIS-I data were cleaned and calibrated using CIAO-4.4. The 
exposure corrected mosaicked image shown in Fig.[TJ(right panel) 
was obtained using the MERGE.OBS script in CIAO. Point sources 
were identified using WAVDETECT using a range of wavelet radii 
between 1 and 16 pixels to ensure all point sources were detected. 

2.3 ROSAT data reduction 

The ROSAT PSPC observations w ere reduced using the Extended 
Source Analysis Software (ESAS, Snowden et a l.lll994l) . and fol- 
lowing the procedure described in lEckert et al.l d2012l) . The result- 
ing exposure corrected mosaicked image for the R37 band (0.4- 
2.0keV) is shown in the right panel of Fig. [TJ When used later to 
derive the density profile the point sources were removed using the 
program DETECT with a constant threshold flux of 0.003 cts/s in 
the R37 band to ensure that the CXB is uniformly resolved across 
the whole field of view. In the mosaicked image in Fig.QJthe point 
sources are left in for display purposes only. 



3 BACKGROUND SUBTRACTION AND MODELLING 

Accurate measurements of the ICM in the cluster outskirts requires 
an accurate knowledge of the components of the X-ray and non X- 
ray background, and an accurate understanding of how these are 
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expected to vary spatially. Suzaku's low earth orbit gives it a low 
and stable non X-ray background which can be reproduced with 3 
percent accuracy using night Earth observations dTawa et al.l2008h . 

3.1 The CXB 

The cosmic X-ray background (CXB) consists of unresolved point 
sources, and is mode l led as a p owerlaw with index 1.4. As de- 
scribed in lBautz et all d2009h and lWalker et alj ( 120121) . our Suzaku 
data alone allow us to remove point sources down to the threshold 
detection limit of Sf^f^lO -13 erg cm" 2 s" 1 in the 2-10 keV 
band, and these were removed during the spectral e xtraction us- 
ing 2. 5 arcmin circular exclusion regions. Following Walk er et al.l 
(2012), we use the two-power-law model for the cumulative flux 
distribution of point sources found using ROSA T, Chandra and 
XMM-Newton observations in iMoretti et al.l j2003l) to find the un- 
resolved CXB flux remaining once the point sources resolved with 
Suzaku have been removed. Taking the total CXB flux in the 2-10 
keV band to be 2.1 8 ±0.13 x 10" 11 erg cm" 2 s" 1 deg -2 (from 
More tti et a j]|2009l using Swift data), and performing the integral, 

Fcxb = 2.18 ± 0.13 x 10" 11 - J {^jg) y-SdS (1) 

we find that following the removal of point sources down to the 
threshold flux Sf" e f afcu the remaining CXB level is 1.87 ± 0.13 
x 10" 11 erg cm -2 s" 1 deg" 2 . 

Outside r = 50 arcmins, our Chandra observations al- 
low point sources to be resolved down to a threshold flux of 
S°J% ndra =1.5x 10~ 14 erg cm" 2 s _1 in the 2-10 keV band across 
the whole ACIS-I field observed. This means that in the outskirts 
the unresolved CXB emission is reduced to 1.37 ± 0.08 x 10" 11 
erg cm -2 s _1 deg -2 . 

To account for the point sources resolved by Chandra in our 
spectral fitting to the Suzaku data, we used the ray-tracing simula- 
tor XISSIM to simulate the emission from the point sources resolved 
with Chandra in the Suzaku pointings for each detector. The spec- 
tra of the point sources extracted from the Chandra data are all in 
good agreement with a powerlaw of index 1 .4, which was used as 
the input spectral model to XISSIM. We simulated exposures 1000 
times the actual exposures to get good statistics. Spectra were then 
extracted from each simulated observation using the annuli shown 
in Fig.Q] and the resolved point source contribution in each annulus 
was modelled using xspec, and included in the background model. 
This means that the only unresolved CXB emission outside r —50 
arcmins is from point sources below 1.5 x 10~ 14 erg cm -2 s _1 . 

Because the CXB consists of unresolved point sources, which 
are not uniformly distributed across the sky, the CXB level deviates 
from the mean when analysing small areas due to varying num- 
bers of unresolved point sources. The expected deviation from the 
average value for a given observed solid a ngle (Q.) resolved to a 
threshold flux S eX ci can be calculated using (Bautz etal. 20fjl), 

^cxb = ^ S ° XCl (^) x S 2 dS (2) 

Our Chandra observations of the outskirts allow the CXB to be 
resolved to a threshold flux ~ 10 times lower than that achiev- 
able with just Suzaku data. Integrating equationf2] we find that this 
reduces the uncertainty in the CXB level due to unresolved point 
sources by a factor of 2.3. 

Due to the large spatial extent of Centaurus, the extraction re- 
gions used in spectral analysis can be made large, which further re- 
duces the expected variations in the CXB level for the regions anal- 



ysed, (as well as significantly reducing the effects of PSF spreading 
between annuli). The expected 1 a variations in the 2-10keV CXB 
flux for the annuli investigated are shown in table [2] Later in sec- 
tion [4/3] we vary the CXB level of the background model through 
these ranges to calculate the systematic error of the CXB level un- 
certainty on the spectral fitting parameters. 

When performing the fits to the cluster emission in the out- 
skirts (outside 50 arcmins) the low temperature of the cluster emis- 
sion meant that the emission above 5 keV was purely from the 
CXB, which allowed us to check that the CXB level was consis- 
tent with the calculated value (see the spectral fits for the outermost 
4 annuli in Fig. lAlb . In all cases the unresolved CXB level obtained 
by letting the CXB norm be a free parameter was highly consistent 
with the value calculated earlier (1.37 ± 0.08 x 10~ n erg cm -2 
s _1 deg -2 for the regions outside 50 arcmins). 

To test the possibility that some of the point sources identified 
with the Chandra observations are actually bright gas clumps, we 
compared the cumulative number counts of point sources (logN- 
logS) with that of the Chandra Deep Field South (CDFS) survey 
field dLehmer et al .120121) . Any excess over the source population in 
the CDFS survey would be evidence for additional bright sources, 
possibly due to gas clumping. The comparison is shown in Fig. 
[5] where we compare the cumulative source number counts in the 
0.5-8. OkeV band observed in our Chandra fields (black points), 
with that observed in the CDFS (blue triangles). Poisson uncertain- 
ties in the total number of sources observed in our Chandra fields 
dominate the uncertainty, and the 1-ct asymmetric confiden ce lim- 
its sho wn were calculated using the formulae presented in lGehrelsl 
dl986h (their equation 9 was used for the upper limits, and their 
equation 14 was used for the lower limits). The source population is 
in complete agreement with the source population from the CDFS, 
and there is no evidence for an excess over the CDFS level which 
would occur if some of the point sources removed were bright gas 
clumps instead of background AGN. 



3.2 Modelling the soft foreground. 

We use the Suzaku region outside 90 arcmin, together with RASS 
background region between 1.5 to 2.5 degree^] (corresponding to 
the region between l.l-1.8r 2 oo). We also use the ROSAT PSPC 
background regions around nearby NGC 4507. All of these back- 
ground regions are well fit by a model of the form phabs*(powerlaw 
+ apec(0.22keV)) + apec(0. 12keV), which consists of an absorbed 
powerlaw to describe the CXB from unresolved point sources, an 
absorbed thermal component at 0.22keV to describe galactic halo 
(GH) emission, and an unabsorbed thermal component at 0. 12keV 
to describe the Local Hot Bubble (LHB). Adding an additional ab- 
sorbed jherm^_ramponent, which is sometimes necessary to fit the 
data dSnowden et alj2008h . does not improve the quality of the fits. 
The temperatures of the galactic halo and local hot bubble were ob- 
tained by spectral fitting to the RASS data, and the metallicities of 
these components were fixed to 1 Zq, while the redshifts of these 
components were fixed to zero. The different background regions 
studied are all found to be reasonably consistent and the best fit pa- 
rameters shown in tablef3] We can use these to gain an understand- 
ing of the expected spatial variation of the soft foreground over the 
cluster, and we later calculate the systematic errors resulting from 



5 The RASS data were obtained from the X-ray background tool at 
http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/xraybg/xraybg.pl 
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Table 2. For each annulus we show the threshold flux for excluding point sources (S excl ), the unresolved CXB flux (F^™ r J so '" ed ), the resolved CXB flux 
(F£, e £°g ved ), and the expected one sigma variations in the unresolved CXB flux each annulus ((^unresolved) m tne 2-10 keV band. 



Annulus 


10-15 


15-20 


20-25 


25-30 


30-35 


35-45 


45-50 


50-60 


60-70 


70-80 


80-90 


90-105 


Sexcl CIO" 14 erg s- 1 cm" 2 deg" 2 ) 
pum-esolved (if)" 12 erg s" 1 cm" 2 deg" 2 ) 

Fcx°B ei (iO -12 er g s" 1 cm ~ 2 deg~ 2 ) 
^unresolved (10~ 12 erg s _1 cm -2 deg -2 ) 


10 

18.7 



3.8 


10 

18.7 



3.3 


10 

18.7 



3.4 


10 

18.7 



3.4 


10 

18.7 



3.3 


10 

18.7 



2.9 


10 

18.7 



3.4 


1.5 
13.7 
2.7 
1.1 


1.5 
13.7 
3.6 
1.1 


1.5 
13.7 
2.4 
1.1 


1.5 
13.7 
7.5 
1.2 


1.5 
13.7 
3.3 
0.8 




S= 0.5-8.0keV flux (erg cm" 2 s ) 



Figure 3. Here we compare the cumulative source number counts obtained 
in our Chandra observations (black points) in the 0.5-8.0keV band with 
the same distribution obtained in the Chandra Deep Field South (CDFS) 
survey (blue points). We find agreement between the two, and there is no 
evidence for an excess in the number counts in our Chandra observations 
which would occur if some of the point sources were bright gas clumps 
rather than background AGN. 



this uncertainty on the background model on the spectral fits in sec- 
tion [43] When performing the background fits the normalisations 
of the LHB, GH and CXB components, and the temperatures of the 
LHB and GH components, were free parameters. 



3.3 Stray Light and PSF spreading 

To calculate the contributions in each annulus from neighbouring 
annuli due to PSF spreading and stray light, we simulated the con- 
tribution in each annulus from all of the other annuli using the 
background subtracted image of the clus ter as input to the ray 
tracing simulator XIS SIM, as was done in IWalker et al.l d2012h and 
Walk er et al.1 d2012ah . The results are tabulated in table [4] where 
each row shows that contribution in each annulus from the annuli 
in the columns. The majority of the emission in each annulus orig- 
inates from that annulus, and most of the external contribution is 
from the annulus immediately inside the annulus in question. 

Due to the roll angle of the observations, stray light from the 



core has been heavily reduced, such that in the outer annulus it 
contributes only 0.3 percent. When performing the spectral fits we 
tested the effect of including light spread from the core and from 
the adjacent annuli, and found this to have no effect on the best fit- 
ting temperatures, densities and metallicities. To perform this test 
we fitted all of the annuli simultaneously and modelled each annu- 
lus as consisting of emission from the annulus itself and all of the 
other annuli. Each annulus was modelled as the sum of the APEC 
component originating from that annulus with those APEC contri- 
butions from the other annuli, weighted by the fractions shown in 
tableg] 



4 SPECTRAL ANALYSIS 

Spectra were extracted in the annuli shown in Fig. [T]between radii 
10-15, 15-20, 20-25, 25-30, 30-35, 35-40, 40-45, 45-50, 50-60, 60- 
70, 70-80 and 80-90 arcmins. These radii were chosen to ensure 
that each annulus contained at least 2000 counts following back- 
ground subtraction. In all cases the annulus width is much larger 
than the PSF of Suzaku (HPD 2 arcmins) to prevent the effects of 
spillage between the annuli. The spectra from each detector (XIS0, 
XIS1, and XIS3) and each editing mode (3x3 and 5x5) were all 
fitted simultaneously for each annulus. 

We first performed projected fits for all of the annuli , mod- 
elling each annulus as an absorbed APEC dSmith et"al] 1200 lh 
component. We fix the column density to the LAB survey 
( Kalb erla et alj 2005t) values for e ach pointing shown in table \t\ 
As in lsimionescu et al.l d201 lh andlUrbanetaLU201ll) . we use the 
abundance tables of Grevesse &Sauvall ( fT998h . 

The projected density profile, obtained from the APEC nor- 
malisations during the spectral fitting, is shown in green in Fig. [4] 
and compared to the projected density profile obtained later in sec- 
tion |4.2| using the ROSAT PSPC data (which have higher azimuthal 
coverage) to ensure the Suzaku results are not biased by looking 
along one strip. We find strong agreement between the Suzaku and 
ROSAT projected density profiles out to 80 arcmins. We find no 
statistically significant cluster emission outside 80 arcmins. 

Deprojected temperature and density profiles were obtained 
by modelling each annulus as the superposition of the ICM from 
the shell corresponding to that annulus with the emission projected 
onto the annulus from the shells exterior to it, with each component 
scaled according to the volume contributing to each annulus. The 
emission from each shell is modelled as an absorbed apec compo- 
nent. This emulates the xspec mixing model PROJCT, but allows 
for the use of background modelling and f or the simultaneous fit- 
ting of data from multiple detectors (as in lHumphrev et al1l2012l ; 
Walke r et al]|2012h . The deprojected densities are derived from the 
APEC normalisations for each annulus obtained during the spectral 
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Table 3. Soft foreground measurements from ROSAT PSPC and Suzaku data. The units of the APEC normalisations are 10- 14 (4tt)- 1 D- 2 (1 + 
z)~ 2 J n e ri[{dV, where Da is the angular size distance (cm), and n e and njj are the electron and hydrogen densities (cm -3 ) respectively, and these 
values are scaled for a circular area of sky of 20' radius (1257 arcmin ). 



Position 0.22 keV APEC norm (GH) 0. 12 keV APEC norm (LHB) 



SUZAKU N offset (90'-105') 


1 7+ ' 5 x 
i./_ 7 X 


10~ 3 


1 5+ 20 


X 


10" 


-3 


RASS 1.8-2.5 degree background 


2 3 +0 ' 4 x 

i-9l :is > 


10" 3 


1 0+ 01 


X 


10" 


-3 


ROSAT NGC 4507 


c 10~ 3 


n -1+0.3 
z - J -0.3 


X 


10" 


-3 



Table 4. Percentage contribution of flux in the rows' annulus from the columns' annulus due to PSF spreading and stray light. 





0'-5' 


5'-10' 


10'-15' 


15'-20' 


20' -25' 


25'-30' 


30'-35' 


35'-45' 


45'-50' 


50'-60' 


60'-70' 


70'-80' 


0'-5' 


94 


5.6 


0.12 


0.039 


0.020 


0.020 


0.020 


0.020 


0.029 


0.029 


0.0099 


0.0099 


5'-10' 


12 


84 


2.2 


0.16 


0.14 


0.12 


0.11 


0.12 


0.13 


0.11 


0.070 


0.029 


10'-15' 


0.61 


7.2 


80 


9.9 


0.34 


0.099 


0.099 


0.14 


0.12 


0.16 


0.13 


0.089 


15'-20' 


0.14 


0.59 


10.0 


80 


7.3 


0.41 


0.22 


0.23 


0.14 


0.14 


0.11 


0.079 


20'-25' 


0.12 


0.19 


0.85 


11.0 


80 


4.7 


0.59 


0.47 


0.38 


0.25 


0.14 


0.12 


25'-30' 


0.11 


0.17 


0.26 


0.55 


6.1 


80 


10 


0.47 


0.29 


0.29 


0.28 


0.20 


30'-35' 


0.31 


0.26 


0.35 


0.37 


0.64 


8.8 


82 


5.4 


0.28 


0.26 


0.19 


0.17 


35'-45' 


0.89 


0.94 


0.94 


0.95 


1.1 


1.6 


7.0 


80 


3.2 


0.66 


0.44 


0.41 


45'-50' 


0.38 


0.34 


0.32 


0.28 


0.25 


0.31 


0.23 


3.4 


87 


6.5 


0.28 


0.26 


50' -60' 


1.5 


1.7 


1.5 


1.6 


1.2 


1.1 


1.2 


0.98 


3.0 


83 


1.2 


0.61 


60' -70' 


0.58 


0.71 


0.81 


0.81 


0.65 


0.78 


0.69 


0.63 


0.59 


4.0 


84 


4.8 


70'-80' 


0.32 


0.58 


0.87 


1.0 


1.0 


1.0 


1.1 


1.0 


0.96 


0.77 


4.9 


86 



fitting, using the volumes of the shells observed in each annulus. 
The spectral fits are shown in Fig. lAll 

The deprojected temperature and density profiles are shown 
in the top two panels of Fig. [5] The temperatures found for the re- 
gion within 0.03-0. 2ri8o with XMM-Newton in iMatsushital hoilh 
are shown by the green points in the top panel of Fig. [5] showing 
excellent agreement with our Suzaku temperatures. The tempera- 
tures f ound for the NW sector using BeppoSAX in Molen di et al.l 
d2002h are shown as the pink points, which again show excellent 
agreement. The red points sho w the temperatures in the NW sector 
obtained using Chandra data in Sanders & Fabian] | |2006|) . 

Using our Suzaku data the metallicity can be constrained out 
to 50 arcmins. The metallicity profile, shown as the black points 
in the bottom panel in Fig. [5] falls from around QAZq in the 10- 
15 arcmin annulus to 0.1 Zq in the 45-50 arcmin annulus. Outside 
of this annulus the metallicity cannot be constrained and it is fixed 
to 0.1 Zq (shown by the blue lines). Therefore, despite the high 
central metallici ty of the Centauru s cluster (reaching 3 Zq in the 
central 30 kpc, iFabian et alj|2005l) the metallicity declines in the 
outskirts and is not higher than metallicity meas urements for other 
cluster s in the outskirts. For the Perseus cluster, Isimionesc u et al.l 
d201 ll) found metallicities arou nd 0.3Zq in t he out skirts, while for 
the less massive Virgo cluster lUrban et"ail d201ll) found a lower 
outskirts metallicity of 0. \Zq. This may suggest that less massive 
clus t ers ha ve lower metallicities in the outskirts. For instance in 
ISunl d2012h it has been found by comparing the abundance profiles 
of samples of galaxy clusters and galaxy groups (their Fig. 4) that 
galaxy groups are iron poorer than galaxy clusters in the region 
0.3-0.7r 50 o. 

The metallic ities found in the region within 0.03-0. 2riso with 
XMM-Newton in iMatsushitd l2oTl|) are shown by the green points 
in the bottom panel of Fig. [5] and agree with the Suzaku metal- 
lities in this region. The metallicities found in the NW sector in 



ISanders & FabiarJ d2006l) are shown as the red points. The pink 
points show the metallicities obtained for the NW sector using Bep- 
poSAX data in Mole ndi et alj |2002h. but we note that th is study 
used the abundance tables of lAnders & Grevessd dl989l) . which 
yield Fe abund ances lower by a factor of a round 1 .4 than the abun- 
dance tables of lGrevesse & Sauvall dl998l) which we use. 

4.1 Entropy and pressure profiles 

The deprojected entropy (K = kT/rit ) is calculated from the 
deprojected temperature and density profiles. The deprojected en- 
tropy profile is shown in Fig. [5] where it is compared to the base- 
line en tropy profile for purely gravitational collapse from lVoit et al.l 
d2005l) , 

K(R)/K 500 = 1.47(r/r 5 oo) 1 - 1 . (3) 

Here we have accounted for the hydrostatic mass estimates being 
biased low by ~ 13 ± 16 percent on average (because these neglect 
the non-thermal pressure support in the ICM), by inc reasing the 
factor in equation [3] from 1 .42 to 1 .47 as described in Pratt et al.l 
d2010h . The baseline profile has b een scaled using the self-similar 
entropy at rsoo (equation[4j as in lPratt et al.l d2010l) . 

K 500 = 106 keV cm' 2 ( —^22 ) ^ ( ±) ^ 

x£( 2 )- 2/3 ft 7 - 4/3 (4) 

using fb=0- 15 as in lPratt et al.1 d2010f) . and using the A/500 value of 
1.2x 10 14 Mq which we find later in our mass analysis in section 

s 

We find that the entropy exceeds the baseline level within 40 
arcmins, and flattens towards it at ~0.5r2oo, wh ich is simil a r to th e 
behaviour found for the REXCESS clusters in lPrattet all feOlfjl) . 
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Outside 50 arcmins (~ rsoo), the entropy is systematically slightly 
below the baseline relation, possibly the results of gas clumping, 
though in the outermost two annuli it is in good agreement with the 
baseline relation. Later, in section [6] we will compare the entropy 
profiles of other clusters explored in the outskirts with Suzaku when 
scaled using self-similar scaling relations. 

The deprojected pressure profile (P = kn e T) is shown in 
Fig. HQ and i s com pared with the universal pressure profile of 
lAmaudetailfcOlOl) . 



P(r) = Ps 

where, 

¥(x) 

(Po, c 5 oo,7, a, /3) 



M soo 



3 x 10 14 Mp 



otj>-\-a'p(x) 



(c 500 Xp[l + (c 5 00^) a ] (/3 - 7)/ ° 
-3/2 



(5) 



(6) 



.403 rift, '\ 1.177, 0.3081, 1.0510, 5.4905) 



q p + a'p(x) = (0.10 + dp) 
x = r/r 500 



1 - 



(x/0.5) 3 



1. + (x/0.5) 3 



(7) 
(8) 
(9) 




after being appropriately scaled by the characteristic pressure, P500, 



Psoo = 1.65 x 10~ 3 h(z) 



8/3 



M500 



-1 2/3 



1 70 1V10 



40 60 
r (arcmin) 



Figure 4. Green shows Suzaku projected densities, black shows ROSAT 
projected densities (for full azimuthal coverage). The position of r200 and 



h 70 keV cm its error is shown by the vertical dashed black and red lines respectively. 



(10) 



The (A/500/3 x 10 14 h 7 " 1 M Q ) Qp+Q p (x) term describes the devia- 
tion from the self similar scaling relation for pressure. From equa- 
tion (8) we see that when x = r/rsoo > 1 the index ap + a' P (x) 
tends to zero, so outside rsoo this term is negligible. 

Whilst the observed pressure profile agrees with the univer- 
sal pressure profile within the spread of the simulations, it is 
systematically slightly higher than the universal pressure profile 
in the outskirts, (as has a lso been found for the Perseus cluster 
in Simionescu et alj|201ll for the fossil group RX Jl 159+5531 
in Humphrey et alj|2012l and for PKS 0745-191 in I Walker etall 



of the entropy profile lying below the baseline entropy profile in 
the outskirts is purely due to the overestimate of the gas density 

<»;..) 



due to clumping. Defining the clumping factor as C 
(Nagai&Lau 201 if), the observed density is overestimated by a 
factor of vC. If the true entropy in the outskirts is given by the 
baseline entropy profile (equation [3j then this means that using the 
observed temperatures and densities at radii R (kT(R) and n e (R)) 
we can determine C(R) using; 

kT(R) 



l2012ah . which further suggests that gas clumping may be having 
an effect on our measured profiles. By causing the gas density to 
be overestimated, gas clumping would cause the gas pressure to be 

overestimated. 

W e can combine the baseline entropy profile of IVoit et al.l 
J2005h and the universal pressure profile of lArnaud et al.l (I20101) . 
appropriately scaled by the self similar entropy (7v"soo) an( l pres- 
sure (P500) respectively, to produce predictions for the temperature 
and density profiles in the outskirts as follows, 

kT(r) = P(r) 2/5 K{r) 3/5 



K(R)/K S00 = 1.47(i?/r 50 o) i 



^50o(ne(7?)/yC(l)) 2 / 3 



(13) 



3/2 



(14) 



n H (r) = (l/1.2)P(r) 3/5 K(ry 



3/5 



(11) 
(12) 



and these predictions are shown as the cyan lines in the tempera- 
ture and density plots in Fig. [5] We see that the temperatures agree 
with this prediction, but the densities are higher than it. This would 
appear to support the suggestion that it is the overestimate of the 
gas density, possibly due to gas clumping, that causes the entropies 
to be slightly lower than the baseline level and the pressures to be 
slightly higher than the universal pressure profile in the outskirts. 

In Fig. [6] we show that if the pressures are corrected by the 
clumping factor needed to make the entropy profile agree with the 
baseline profile in the outskirts (shown later in Fig. Ill) , then they 
agree better with the universal pressure profile. These clumping 
factors have been calculated using the assumption that the cause 



x/CCR) = n e (R) (W^soo) 11 J|gy) 

If gas clumping is present and the gas clumps are in pres- 
sure equilibrium with the surrounding ICM, then they would be 
expected to be colder than the surrounding IC M, and therefore bia s 
the temperature low (this is also discussed in lUrban et al. 1 1201 lh . 
However we have found that for Centaurus the temperature profile 
in the outskirts is in good agreement with the temperature profile 
obtained by assuming the universal pressure profile and the baseline 
entropy profile, and that the deviations from the universal pressure 
profile and the baseline entropy profile appear to be purely the re- 
sult of a density excess. This may indicate that the gas clumps are 
not in pressure equilibrium with the surrounding ICM. For instance 
it is likely that any gas clumps will be falling into the cluster and 
moving through the ICM, so ram pressure will help to support the 
clumps, and this may be the dominant support mechanism. 



4.2 Comparing to ROSAT data 

Following the approach o f lEckertetal] | |2012h we convert the sur- 
face brightness profile from the ROSAT PSPC mosaic into a pro- 
jected density profile. As described in section 12.31 point sources 
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40 60 
R (orcmin) 



Figure 5. First panel: Deproj ected tempera ture profile. The green points 
are the temperatures from Matsushita fcpl lb obtained with XMM-Newton, 
the pink points are the temperatures obtained for the NW sector in 
Molen di et ailfeOOA using BeppoSAX, and the red points are the tempera - 
tures obtained in the NW sector with Chandra in Sanders & Fabian (2006). 
The cyan lines shows the temperatur es and densities predicted by com- 
bining the baselin e entropy profile of IVoit et aD J2005h and the universal 
pressure profile of Arn aud et al Second panel: Deprojected density 

profile, Third panel: Deprojected entropy profile. The dashed line shows 
the r powerlaw relatio n predicted for p urely gravitational hierarchical 
structure formation (from IVoit et all d2005l) ). scaled using the self similar 
entropy, if 500 from equ ation [4] Fourth panel: Metallicity profile. Green 
points are from Matsushita ( 201 lj) usi ng XMM-Newton; pin k points are for 
the NW sector using BeppoSAX from Mo lendi et al.l ^2 002); the red points 
are for the NW using Chandra from Sanders & Fabian (2006) . The blue 
lines show where the metallicity has been fixed in the outskirts to 0. 1 Zq . 
In all of the panels the solid black lines show the systematic errors calcu- 
lated in section 1431 The dashed vertical black line shows the value of r2oo 
calculated in section [5] with the vertical dashed red lines showing the error 
inr 2 oo- 




Figure 6. Deprojected press ure profile, which is consistent with the univer- 
sal pressure profile found in lArnaud et al.l (12010'). shown as the solid black 
line inside rsoo (where it is based on observations), and as a dotted line 
outside rsoo (where it is based on simulations). Th e uncerta inties on the 
universal pressure profile in the outskirts shown in lLapi et al] j2012l) are 
shown by the cyan lines. The red points show the effect of correcting the 
pressures by the clumping values needed to make the entropy profile agree 
with the baseline profile in the outskirts. This clumping correction lowers 
the pressures and brings them into even better agreement with the universal 
pressure profile. 



were removed using the program DETECT with a constant thresh- 
old count rate of 0.003 cts/s in the R37 band to ensure that the 
CXB is resolved u niformly across the whole field of view, as in 
lEckert et all J2012h . The background regions of the offset pointing 
of NGC 4507 shown in Fig.[T]were used to obtain the background 
level. The projected density profile was obtained using the ROSAT 
PSPC response and the Suzaku temperatures for each annulus to 
convert the count rates in each annulus into an apec normalisation. 

As the Suzaku observations have small azimuthal coverage, it 
is important to understand whether they are representative of the 
cluster as a whole. We find that the projected density profiles are 
highly consistent, as shown in Fig. [4] 



4.3 Systematic errors 

We need to calculate the systematic errors on the deprojected tem- 
perature, density and entropy profiles resulting from the uncertainty 
of the X-ray background parameters (the GH, LHB and CXB nor- 
malisations) which we have quantified in section ED and the un- 
certainty in the NXB level, which from iTawa et al.l ((2008) is ± 3 
percent. To achieve this we produce 10000 reali sations of the back- 
ground mod el and the N XB level (as done in Walke r etal.ll2012l 
and lWalker et al.|[20 12ah . allowing all of the background parame- 
ters to vary simultaneously within their variances calculated earlier 
(and taking into account the covariance between background model 
parameters), and perform the deprojection in xspec for each real- 
isation. This allows the uncertainty on the background model to 
be folded through the deprojection, resulting in a complete prop- 
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agation of the errors. This allows a more realistic estimate of the 
systematic error to be obtained than is achieved by varying only 
one background parameter and leaving the others fixed to their best 
fit values. 

From the distribution of the 10000 profiles we can find the 1 
a systematic errors, and these are shown by the solid lines in Fig. 
[5] These systematic errors also include the effect of varying the 
contamination on the optical blocking filter by ±10 percent, which 
was estimat ed by using the ftool X ISCONTAMICALC to modify the 
ARFs (as in lAkamatsu et alj|201 11) . and was found to be a negligi- 
ble effect much smaller than the statistical errors. The systematic 
errors also include the effect of varying the column density by ± 
20 percent (the range of the values from the LAB survey for the 
observations as shown in table [T), and of varying the metallicity in 
the range 0.0-0. 3Z in the outer three annuli when performing the 
fits. In all cases the total systematic errors are less than or equal to 
the statistical errors. 



dNavarro et al.lll997h . which is well motivated by numerical sim- 
ulations of gravitational collapse, 



p(r 



Po 



(16) 



r/r s (1 + (r/r 8 )) 
po = 200p c cloo/3(ln(l + C2oo)-C2oo/(l + C2oo)) (17) 

W e follow the same mass analysi s method as Walker et al.l 
which is based on that used in ISchmidt & Ailed J20071) . 
We use the deprojected density profiles to predict the temperature 
profile assuming hydrostatic equilibrium and the NFW profile for 
the total mass, moving inwards from the outermost annulus (the re- 
sults are unchanged if the method starts from the innermost annulus 
and moves outwards). We use the XMM-Newton temperatures near 
the core to better resolve the central regions. 

The best fitting mass profile produces the best fitting temper- 
ature profile which minimises the % 2 statistic 



4.4 Large scale motions 

To search for the possibility of large scale gas motions, we di- 
vided the ROSAT image b y the azimuthal averag e surface bright- 
ness profile (as done in ISimionescu e t al. 2012 for the Perseus 
cluster). We find no significant features in the outer regions. In 
the central regions, the only significan t features we find are those 
which have already been studied in IChurazov et"all d 19991) (their 
Fig. 3, right panel), an d which our direction of study has avoided. 
IChurazovetai1 ( ll999l) found an enhancement around 20 arcmins to 
the south east of the core, and that this enhancement corresponds to 
one of the subgroups in Centaurus (Cen 45, which contains the sec- 
ond brightnest galaxy, NGC 4709). The ASCA temperature map 
of the central regions also indicated that the enhancement in sur- 
face brightness coincides with an enh ancement in temperature to 
the south east. IChurazov et al.l d 19991) concluded that this higher 
temperature arose due to the heating caused by the interaction be- 
tween the Cen 45 subgroup as it merges with the main cluster. This 
merging activity may be responsible for the c old front to the wes t 
of the core studied in depth with Chandra in iFabian et al.l l l2005h . 

by p roducing a slosh i ng mo tion o f the core. 

iTakahashi etal] J2009I) and lLovisari et al. | ([201 ]]) produced 
temperature maps of the central regions of Centaurus with full 
azimuthal coverage using XMM-Newton data and found that the 
north western sector has a slightly lower temperature than the other 
directions, however these observations only extend out to 12 ar- 
cmins (160 kpc=0. 14r2oo). This was a lso found in the Chandra 
temperature map of Lagana et al. d2010l) . but this only extends out 
to 50 kpc (0.05r2oo). Since the surface brightness is only weakly 
dependent on the temperature, any azimuthal variations in temper- 
ature will have a negligible effect on the observed gas mass. 



5 MASS ANALYSIS 

Assuming the ICM to be in hydrostatic equilibrium, and assuming 
the cluster i s spherically symmetr ic, the total mass within radius, r, 
is given bv dVildilininet1ifll2006h 

M« r) = -3.68 x lO 13 M T(r)r ( + ^) (15) 

\ a In r a In r J 

We assume the total mass is described by an NFW profile 



2 \ (7/calculatcd,z ^actual,?) , * n\ 

x =2- 4 (18) 

i 1 actual,! 

As described in Allen et al 1 (E008h . because this method does 
not involve the use of parametric fitting functions for the tempera- 
ture and gas density profiles, it avoids the strong priors that these 
place on the mass determination and which complicate the interpre- 
tation of results (and can lead to an underestimate of uncertainties). 
We stress that there is no robust prior expectation for the form the 
temperature and density profiles in the outskirts of clusters, and so 
the non-parametric approach we employ is important to ensure the 
results are not severely restricted by priors. 

To propagate the uncertainty of the density profile into the 
measure ment of the gas mass fraction and the NFW parameters, 
we follow lWalker et all d2012ah and repeat the mass analysis 10000 
times using different realisations of the density profile distributed 
by the combined statistical and systematic errors shown in Fig. 
[5] For each iteration the cumulative gas mass fraction is calcu- 
lated and the NFW best fitting parameters. The best fit values, to- 
gether with the 1 a errors were then calculated from the result- 
ing distributions. We find that the NFW profile describes the total 
mass profile well. We find C2oo=5.9l*'|, r s =191 Igjkpc, r2oo = 
1130t54 5 kpc= 84ttoarcmins and M 20 o = 1.6±g;f x 1O 14 M , 
and the fitting statistic is x 2 /d.o.f = 12.5/12. 

To co mpare the calculated mass value with the M-T scaling 
relation of Arna ud et al.l | |2005|) . we found the spectroscopic tem- 
perature in the region 0.1 < r < 0.5r2oo, which is fcT(0.1 < 
r < 0.5r 20 o) = 3.0lg;i keV > ( these 

errors include the system- 
atic errors in the background modelling). Using the scaling relation 
M 2O o/lO 14 M = 5.34 ± 0.22 x (fcT/5keV) 1 ' 72±0 - 10 //i(2) we 
find M^ naud =2.2tol x 1O 14 M , which is slightly higher but in 
reasonable agreement with our calculated mass giv en the scatter 
around the best fit relation for the M200-T relation in lArnaud et al.l 
<2005h . 

Because the clusters studied in lArnaud et al] (2005) were 
mostly only studied out to a maximum radius of ~0.6 r2oo, the 
M200 values used to obtain the M200 — T scaling relation were 
derived by extra polating out the bes t fitting NFW models to r2oo • 
As described in lArnaud et al.l d2005l) . the quality of the powerlaw 
best fit decreases during this extrapolation, with the null hypothesis 
probability for the whole sample decreasing from 0.32 at 5=2500 
to 0. 07 at (5=200. In Fig.[7J we plot Centaurus on the M 20 o - T plot 
from lArnaud et al.l d2005l) as the red point, and we see that it is in 
reasonable agreement given the scatter around the best fit relation. 
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Figure 7. The Mian - T plot from lArnaud et al. (2005), with our value for 
the Centaurus cluster ad ded as the red poin t. The solid line shows the best 
fit powerlaw relation in Arnaud et al. (2005j) fo r all of the clusters, w hile the 
dashed line is the best fit powerlaw relation in lArnaudetai1d2005l) for the 
clusters with fcT(0.1 < r < 0.5r 2 oo) > 3.5keV. 



The gas mass fraction profile is shown in Fig. [8] and is found 
to rise to the mean cosmic baryon fraction at the highest radius 
observed. As shown in Fig. [8] the gas mass fraction in the out- 
skirts agrees with that found using Planck data for a st acked sam- 
ple of 62 clusters in [pianck Collabor ation et all d2012l) (when the 
IVikhlininet all 2006 temperature profile is assumed when obtain- 
ing the density profiles from the Planck data). In Fig. |8]we show 
the effect of correcting the gas mass fraction by the clumping fac- 
tors calculated for Centaurus in section [6] which are required to 
make the entropy profile agree with the baseline entropy profile in 
the outskirts. We see that performing such a correction reduces the 
gas mass fraction so that it agrees with the expected hot gas frac- 
tion (which takes into account that 12 percent of the baryons should 
be in stars). In Fig|9]we see that the clum ping corrected gas m ass 
fraction agrees well with the predictions of lYoung et al.l (|20 1 lh for 
clusters in the temperature range 2.5-5.0 keV. 



6 UNIVERSAL ENTROPY PROFILE 

In I Walker et al.l d2012bl) we reported that the entropy profiles of 
clusters explored to r2oo with Suzaku, XMM-Newton and Chan- 
dra have the same shape, flatten i ng off from a powerlaw relation 
above ~0.5r2oo- In lWalkere t al. (2012b) the entropy profiles were 
scaled by their value at 0.3r2oo, where all of the profiles demon- 
strated a powerlaw increase. Greater insight can however be gained 
by scaling the entropy profiles by the self-similar entropy at rsoo 
(K500), as this allows both the shape and normalisation of the en- 
tropy profiles to be compared with theoretical expectations. Similar 
flattening has for instance been observed i n the REXCESS sample 
of clusters obtained with XMM-Newton dPratt et alfoiol) . how- 
ever the entropy is enhanced within rsoo, and the flattening only 
acts to bring th e entropy back in to agreement with the baseline en- 
tropy profile of IVoit etaL at rsoo- 

To scale th e entropy profiles we follow the approach of 
IPratt et all fa) 10). and scale the entropy profiles by the self-similar 
entropy at 7-500 (equation |4j, and we compare this to the base- 
line entropy profile assuming only gravitational physics obtained 
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Figure 8. The cumulative gas mass fraction profile is shown as the solid 
black line with its error shown by the solid cyan lines. The dashed black 
and cyan lines show the effect of performing a clumping correction to 
make the entropy profile in the outskirts agree with the baseline profile. 
The horizontal dashed black line sh ows the mean cosmic baryon fraction 
of 0.167 found using WMAP data in Kom atsu etafl feplll ). while the hor- 
izontal dashed green line shows the expected hot gas fraction when we ac- 
count for 12 percent of baryons being in stars. The vertical dashed black 
line shows r200, and its error is in red. The blue lines are the upper and 
lower limits of the gas mass frac tion outside 0.6r200 found using Planck in 
IPlanck Collaboration et al. 1 20l2) for a stacked sample of 62 clusters (using 
the Vikhlinin et al. 2006 temperature profile), which agrees with the Cen- 
taurus profile in the outskirts. 
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Figure 9. The same as Fig.[8]but here we compare the clumping corrected 
gas mass fraction profile (dashed black line with erro rs in cyan) with the 
expected range in f ga s obtained in the simulations of Youn g~etal] feplll) 
for clusters with an average T in the 2.5-5.0keV range, shown by the green 
lines. 



in l Voit et"al] d2005h (equation[5]l. The clusters used are tabulated in 
table [5] The M500 values are taken from the values q uoted in the 
paper s listed in table \5\ however for the Virgo cluster dUrban et al.l 
I201 lh no mass analysis had been performed, so M500 is calculated 
using the scaling relation o fl Arnaud et al.l d2005h a nd the mean tem- 
perature of 2.3keV reported in lUrban et al.l(l201ll) . 



© 0000 RAS, MNRAS 000, 000-000 



12 S. A. Walker et al. 



Table 5. Sample o f galaxy cluster outskirts observations used. Masses marked with an asterisk were calculated using the M500 — T scaling relation of 
lArnaud et alj d2005h - Masses marked with t were taken from the values measured in lArnaud et al.|j2005l) for those clusters. 
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Figure 10. Entropy profiles of clusters explored with Suzaku, XMM- 
Newton and Chandra in the outskirts, scaled by the entropy at r^oo pre- 
dicted by self-similar scaling r elation s . The solid green line shows the base- 
line entropy profile from IVoit et alj 12005) (equation [3} calculated using 
only gravitational physics. The black line shows the me dian entropy profile 
from the REXCESS cluster sample in lPratt et alJbOlCh . 



In Fig. [TO] we find that the entropy within rsoo = O.659r2oc0 
is in excess of the baseline level for most of the clusters, and in 
general agrees with the baseline level at rsoo- The be haviour in- 
side r sop roughly agrees with the median profile from IPratt et all 
(2010), shown by the solid black line, though we note that the 
REXCESS sample contained many lower mass systems which have 
acted to increase the median entropy level compared to our sample 
of mostly massive clusters. The low mass group RX Jl 159+5531 
and the Virgo cluster demonstrate the largest entropy excess. We 
note that the Virgo entropies are in reasonable agreement with the 



° The relation rsoo = 0.659r200 assumes an NFW profile and the m ean 
concen tration parameter from the sample of Pointecouteau et~al] i2005l) . as 
used in lPratt et alj <20ld) 



Figure 11. Profile of the factors by which the gas density would need to 
be overestimated if gas clumping is the sole cause of the measured entropy 
profiles lying below the baseline entropy profile in the outskirts. 



baseline entropy profile in the outskirts (though the scatter for the 
Virgo datapoints is large). 

However outside 7-500 the entropies are systematically below 
the baseline prediction using only gravitational physics. One pos- 
sibility is that gas clumping (which causes the gas density to be 
overestimated) is causing the entropies outside rsoo to lie below 
the baseline entropy profile. In Fig.fJTJwe show the factor by which 
the gas density needs to be overestimated in order for the entropy 
to agree with the baselin e profile using equati onll4lfor each cluster. 

The simulations of iNagai & Lad ( 1201 II) found the predicted 
clumping level to be mildly mass dependent, with the more mas- 
sive clusters having slightly higher clumping factors (thus causing 
a greater underestimate of the entropy). To search for signs of this 
we plot in Fig. [T2] the mass dependence of the entropy decrement 
below the baseline level (Kvoit / Kobserved) near r2oo for the clus- 
ters studied. Due to the large errors and scatter, and the small sam- 
ple size, there is no statistically significant correlation between the 
entropy decrement below the baseline level and the cluster mass. 
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Figure 12. Plotting the decrement of the entropy below the Voit baseline 
entropy profile at T2oo against cluster mass. There is no statistically signif- 
icant correlation. 



Further insight can be obtained by comparing the scaled tem- 
peratures and densities measured in the outskirts (outside 0.5r2oo) 
with the temperature and density profiles that result when both 
the baseline e ntro py profile and the universal pressure profile of 
Arnaud et al. N2O10l) are assumed to hold in the outskirts. Scaling 
by the self similar entropy and pressure we can write the expected 
temperature profile as, 

(P(r)/P 500 ) 2/5 (K(r)/K 500 ) 3/5 



kT (r) 



(19) 



and the expected hydrogen density profile as, 



scaled 



(r) = (l/1.2)(P(r)/P 50 o) 3/5 (^(r)/ J rr 50 o)- 3/5 (20) 

where K(r) / K500 is given by equation [3] and P(r)/Psoo is given 
by equation|5] Outside r 50 o the (M 50 o/3 x 10 14 hf 1 M Q ) ap+ct p (x) 
term in the pressure profile is negligible. 



We scale the observed temperatures by P 500 



2/5 



-3/5, 
500 



-3/5 
500' 



^500 an( l me 



and compare them to equa- 
tions 1 191 and [2Uj in tig. [T3] No density pr ofiles are p r esente d in 
iBonamente etaU ( l2013l) (for Abell 1835) or lBautz et afl l l2009h (for 
Abell 1795), so these cannot be shown. 

We see that the temperatures generally agree with the pre- 
dicted temperature profile out to r2oo, but the temperatures outside 
this are lower. These low temperatures are the dominant cause of 
the entropy decrement compared to the baseline profile for the 4 
clusters in question (A1689, A2142, A1835 and PKS 0745-191). 
The good agreement with the predicted temperatures within r2oo 
seems to indicate that temperature biases due to cold gas clumps 
are not significant within r2oo, which can be explained if the gas 
clumps are not in thermal pressure equilibrium with the surround- 
ing ICM but are mostly confined by ram pressure as they move 
through the ICM. Outside r2oo it is possible that the decrease in 
temperature below the predicted level is at least partly the result of 
cold gas clumps biasing the temperature low. 

Most of the measured densities lie above the predicted density, 
for which gas clumping may be the cause. This overdensity appears 
to be the cause of the flattening of the entropy profile for Perseus, 
Hydra A and Abell 2029, which all have temperatures in agreement 
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Figure 13. Comparing the self similar scaled temperatures (top panel) and 
densities (bottom panel) outside 0.5r2oo wi th the predictions obtained by 
assuming that the baseli ne entropy profile o f Voit et al. (2005) and the uni- 
versal pressure profile of Arnaud et al. (20 13 both apply in the outskirts. 



with the predicted temperature (i.e. for these clusters the entropy 
decrement does not occur because the temperatures are too low, 
but because the densities are too high compared to predictions). 
By contrast, for A2142 and A1689, the densities are in reasonable 
agreement with predictions in the outskirts, and the entropy decre- 
ment in the outskirts of these clusters appears to be caused by the 
temperatures being too low compared to predictions. 
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6.1 Comparison with Eckert et al. 2013 

Recently, lEckert et al.l d2013al) combined Planck pressure profiles 
and ROS AT PSPC density profiles to derive the entropy profiles for 
a sample of clusters and claimed that the entropy profiles of cool 
core clusters agree with the baseline entropy profile outside r2oo- 
Their results for cool core (CC) and non cool core clusters (NCC) 
are shown overplotted on the Suzaku results as the blue and red 
shaded regions respectively in Fig. [14] A complete investigation 
into the discrepan cy between the Suza ku results outside r2oo with 
those presented in Ec kert et al, (2013a) is beyond the scope of this 
paper, however here we briefly disc uss some pos s ible cau ses. 

To obtain the entropy profiles, lEckert et al.l d2013al) fitted the 
Planck pressure and ROS AT density profiles with functional forms, 
which means that a functional form was assumed for the entropy 
profile. The priors placed on the entropy profile by the assump- 
tion of this functional form, and the degrees of modelling freedom 
available, were not fully demonstr ated. 

As can be seen in Fig C.2 in lEckert et al. (2013a). the uncer- 
tainties on the raw ROS AT density profiles increase dramatically 
outside rsoo, and the functional forms which are fit to these pro- 
files appear unable to fully explore the errors in the outskirts. This 
leads to the error envelope of the functional forms (green shaded 
regions in Fig C.2) significantly underestimating the true density 
errors outside rsoo • It is unclear how sensitive the fits to the ROS AT 
density profiles are to the data points outside rsoo, and no statistic is 
presented to indicate how sensitive the fits are to the outer regions. 
It is therefore unclear to what extent the functional form fitting is 
controlled by the central regions where the data quality is higher. 

The underestimate of the errors through th e use of functional 
forms is evident in Fig. 3 in Eck ert et"aH d2013ah where the temper- 
ature profiles derived from the functional form fitting (solid green 
regions) are compared to those derived using non-parametric de- 
projection (red triangles). We see that outside r2oo the use of the 
functional form causes the temperature errors to be underestimated 
by at least a factor of 3. This is also evident in the gas mass frac- 
tion profiles derived from the same data in Fig. 2 (left panel) of 
lEckert etal.N2013bh . where the use of functional forms underes- 
timates the errors outside r2oo compared to the non-parametric 
method by at least a factor of 3. 

In addition, unlike in iPratt et alj fcoid) where the individual 
entropy profiles for each cluster are displayed to give an indication 
of the scatter around the mean REXCESS entropy profile, the scat- 
ter around th e mean entropy profi le for each individual cluster is 
not shown in lEckert et al.1 d2013al) . so it is unclear what the range 
of the entropy profiles is. There is for instance a large range in the 
gas mass fra c tion pro files for individual clusters shown in Fig. 1 of 
lEckert et al.l d2013bl) . which is much larger than the mean profile 
(green shaded area ) they claim f r om th eir functional form fitting. 
The simulations of iNagai & Laul fcOlfh found that the amount of 
gas clumping can vary significantly between individual clusters, re- 
sulting in a large scatter around the median entropy profile in the 
outskirts. 



7 SUMMARY 

We have explored the thermodynamic properties of the ICM of the 
outskirts of the Centaurus cluster, reaching out to 95 percent of the 
T2oo value we measure along a strip to the north west. The den- 
sity profile we have found with Suzaku agrees well with the den- 
sities obtained using ROSAT PSPC observations which sample the 
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Figure 14. Same a s Fig. 1101 but with the entropy profiles obtained in 
lEckert et alj d2013al) overplotted. Blue i s for CC clusters an d red is for NCC 
clusters (the same colour notation as in Ec kert et al .l2013al) . 



whole azimuth, indicating that our densities are not biased by look- 
ing along only one strip. 

The entropy profile demonstrat es the same centr al excess re- 
ported for the REXCESS clusters in Pratt et al. (2010), and agrees 
with the baseline entropy profile of IVoit et al.l (120051) in the out- 
skirts, but is systematically slightly below it (Fig. |5}. We find that 
the pressure profile agrees with the universal pressure profile of 
lArnaudetaf] ( l2010l) . but lies systematically slightly above it in the 
outskirts. The lower entropies and higher pressures in the outskirts 
may be the results of gas clumping in the outskirts causing the gas 
density to be overestimated. This conclusion is supported by the 
fact that while the temperatures in the outskirts agree with the tem- 
peratures obtained by combining the baseline entropy profile and 
the universal pressure profile, the densities lie above this prediction 
(Fig.[5}. 

The gas mass fraction profile we measure rises to the mean 
cosmic baryon fraction near r2oo, and agrees wi t h the values 
obtained with Planck in IPlanck Collaboration et al.1 d2012l) at the 
largest radius studied. Correcting for the gas clumping needed to 
make the entropy profile agree perfectly with the baseline profile in 
the outskirts causes the gas mass fraction to agree with the expected 
hot gas fraction (which accounts for 12 percent of the baryons being 
in stars), and also brings th e gas mass fraction profile into agree- 
ment with the predictions of You ng et al. I d201ll) . 

(Fig- ID- Perform- 
ing this clumping correction also brings the pressur e profile into 
better agreement with the universal pressure profile of lArnaud et al.1 
d201Ch (Fig.©. 

We find that the total mass profile is well described by an NFW 
profile with c 2 oo=5.9tJ;®, r s =19ll™kpc, r 20 o = 1130l^' 5 kpc= 
84t4;oarcmins and M 20 o = 1.6±o f x I0 14 M©. The derived mass 
is in reasonab le agreement with the M200 — T relation found in 
lArnaud et al.l (2005) given the scatter around the best fit relation, 
and lies below the relation, in agreement with the steepening of the 
M200 — T relation away from the self-similar Al cc J 13 / 2 relation 
when lower mass clusters are included. 

We have furthered the analysis of the collective properties 
of the entropy profiles of clusters explored in the outskirts with 
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Suzaku originally reported in Wal ker etal.ld2012bl) . When the en- 
tropy profiles are scaled by the self-similar entropy, 7\"soo, we find 
that the clusters g enerally demonst rate an excess above the baseline 
entropy profile of Voit et al. (200 5b within rsoo, as has been found 
for the RE XCESS clusters i n IPratt et alj d2010h and for a sample 
of groups in I Sun etHll d2009h - The entropy profiles then flatten and 
agree with the baseline entropy profile at around rsoo. However 
outside rsoo the entropy profiles tend to lie below the baseline en- 
tropy profile, indicating that non-gravitational processes are present 
outside rsoo which lower the entropy. One possibility is that gas 
clumping is respo nsible, (as has been shown in the simulations of 
Naga i & Laull201 ]]), and the overdensities required to bring the en- 
tropies into agreement with the baseline entropy profile are calcu- 
lated in Fig. II II 

We have compared the scaled temperature and density pro- 
files in the outskirts of these clusters (outside 0.5r2oo) with the 
profiles obtained by assuming both the baseline entro py profile of 
IVoit et al.1 {2005) and the universal pressure profile of lArnaud et al.1 
(2010) (Fig.|13t. The temperatures agree with this prediction in the 
range 0.5-1.0r2oo, however the 4 measurements which have been 
obtained outside r 20 o (for PKS 0745-191, Abell 1835, Abell 2142 
and Abell 1689) all lie below the predicted temperature, causing 
the entropy to be lower than the baseline level. The density profiles 
tend to lie above the predicted level outside 0.6r2oo. possibly as a 
result of gas clumping, which acts to reduce the measured entropy. 
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APPENDIX A: 



The outskirts of Centaurus 17 
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Figure Al. Spectral fitting for the cluster emission. In reading order we show the spectra from the annuli between 10^- 15' , 15' -20', 20 / -25 / , 25 / -30 / , 30' -35', 
35 / -45 / , 45 / -50 / , 50' -60', 60 / -70 / and 70' -80'. The data from each detector and each editing mode were fitted simultaneously and are added here for display 
purposes only. The black lines through the points represent the best fits (background plus cluster emission), while the lower black line shows the cluster 
emission. The green line shows the X-ray background level. 
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